The elusive Heisenberg limit in quantum enhanced metrology 
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Quantum precision enhancement is of fundamental importance for the development of advanced 
metrological optical experiments such as gravitational wave detection and frequency calibration 
with atomic clocks. Precision in these experiments is strongly limited by the 1/yN shot noise 
factor with N being the number of probes (photons, atoms) employed in the experiment. Quantum 
theory provides tools to overcome the bound by using entangled probes. While in an idealized 
scenario this gives rise to the Heisenberg scaling of precision 1/N, we show that when decoherence is 
taken into account, the maximal possible quantum enhancement amounts generically to a constant 
factor rather than quadratic improvement. We provide efficient and intuitive tools for deriving the 
bounds based on the geometry of quantum channels and semi-definite programming. We apply these 
tools to derive bounds for models of decoherence relevant for metrological applications including: 
dephasing, depolarization, spontaneous emission and photon loss. 
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Quantum enhanced metrology aims to exploit quan- 
tum features of atoms and light such as entanglement, 
for measuring physical quantities with precision going 
beyond the classical limit 1 . A prominent example is 
that of an optical interferometer, where interference of 
photons at the output port carries information on the 
relative optical path difference between the interferome- 
ter arms. When standard laser light is used, the observed 
results are compatible with the claim that "each photon 
interferes only with itself" 5 and the whole process may 
be regarded as sensing with N independent probes (pho- 
tons) . Parameter estimation with N independent probes 
yields the l/y/~N standard scaling (SS) of precision 6 . En- 
tangling the probes however, can in principle offer a 
quadratic enhancement in precision, i.e. the 1/N or 
Heisenberg scaling (HS) 7-11 . Such strategies have been 
experimentally realized in optical interfcrometry 12-16 
with exciting applications in the quest for the first di- 
rect detection of gravitational waves 17,18 . Moreover, the 
same quantum enhancement principle can be utilized in 
atomic spectroscopy 19,20 where the spin-squeezed states 
have been employed for improving frequency calibra- 
tion precision 21-24 . Alternative approaches to beat the 
SS without resorting to quantum entanglement include 
multiple-pass 25,26 and non- linear metrology 27,28 . 

Unfortunately, both the theory 29-33 and experiments 34 
confirmed the fragility of the above schemes when noise 
sources such as decoherence are considered, and it has 
been rigorously shown for particular models that asymp- 
totically with respect to TV, even infinitesimally small 
noise turns HS into SS, so that the quantum gain 
amounts to a constant factor improvement 3 . 

In this paper we develop two methods which allow to 
extend these partial results to a broad class of decoher- 
ence models, and in particular to obtain fundamental 
bounds on quantum enhancement for the most relevant 
models encountered in the quantum metrology literature. 
The two methods complement each other in terms of pro- 
vided intuition and power. 
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Figure 1 J Quantum metrology and the classical simulation 
idea, a, General scheme for quantum enhanced metrology. N- 
probe quantum state fed into N parallel channels is sensing an 
unknown channel parameter tp. An estimator tp is inferred from 
a measurement result on the output state, b, Classical simula- 
tion of a quantum channel. The channel is interpreted as 
a mixture of other channels Ax, where the dependence on tp is 
moved into the mixing probabilities p v (X). 



The first one elaborates on the idea of classical simu- 
lation (CS) 38 and provides a bound based solely on the 
geometry of the space of quantum channels (SQC). When 
applicable, it gives an excellent intuition as to why the 
HS is lost in the presence of decoherence, but fails to yield 
useful bounds for some relevant decoherence models. 

The second one is based on the channel extension (CE) 
method 39 and requires the optimization over different 
Kraus representation of a quantum channel. However, 
unlike the earlier work 3 which involved making an edu- 
cated guess about the appropriate class of Kraus repre- 
sentations, our bound can be cast into an explicit semi- 
definite optimization problem which is easy to solve even 
for complex models. The power of the method is demon- 
strated by obtaining new bounds for depolarization and 
spontaneous emission models, and re-deriving asymptotic 
bounds for dephasing and lossy interferometer with un- 
precedented simplicity. 



BOUNDS ON PRECISION IN QUANTUM 
ENHANCED METROLOGY 

The typical quantum metrology scenario is as follows 
(see Fig. la). An ensemble of N quantum systems un- 
dergo in parallel the same transformation A v , which de- 
pends on an unknown physical parameter tp. The output 
state is measured and the outcome is used to compute 
an estimate Cp of the parameter ip as summarised below 

ip^A® N [p N ] —tip. (1) 

The task is to find the optimal (possibly highly entan- 
gled) input state p N and the most effective measurement 
strategy in order to minimize the estimation error Ay>jv- 
Note that since the decoherence process is assumed to act 
independently on each of the probes, the global channel 
is described by the tensor product A® N . 

We pursue the estimation problem by restricting our 
attention to estimators which are unbiased in the neigh- 
bourhood of some fixed parameter value, for which the 
quantum Cramer-Rao bound holds 

A<p N > , 1 =, (2) 
y/F Q [A$ N \p»]] 

where Fq is the quantum Fisher information (QFI) 9 . 
Maximization of QFI over input states p N sets the limit 
on the achievable quantum enhanced precision: 

A m > - * F N = max Fq [A® n [p N ] ] . (3) 

VJ-N P N 

The states that maximize the QFI and yield the 
HS for decoherence-free unitary channels are typically 
highly entangled: the GHZ state in the case of atomic 
spectroscopy 40 , and the NOON state in the case of optical 
interferometry 41 . In the presence of decoherence, the op- 
timal input states do not have an intuitive form and the 
maximization of the QFI (with rising N) becomes hard 
even numerically 20 ' 29,30 . 

The typical behaviour of the estimation uncertainty 
in the presence of decoherence is depicted in the log- 
log scale in Fig. 2, showing that asymptotically in N the 
quantum gain amounts to a constant factor improvement 
over the standard 1/\N scaling achievable with indepen- 
dent probes. The key result of this paper is to provide a 
method for a general and simple calculation of this con- 
stant factor improvement. 

CLASSICAL SIMULATION 

To understand the idea of CS, we need to think of 
quantum channels in a geometrical way 42 . A quantum 
channel is a completely positive, trace preserving map 
acting on density matrices. The space of all such trans- 
formations is convex: if A, A' are two channels, then 
pA + (1 — p)A' can be realised by randomly applying 
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Figure 2 Estimation precision in presence of decoherence. 

Log-log plot of a generic dependence of quantum enhanced pa- 
rameter estimation uncertainty in the presence of decoherence as 
a function of the number of probes used. While for small number 
of probes the curve for achievable precision follows the Heisen- 
berg scaling, it asymptotically flattens to approach const/V^V 
dependence. The "const" represents the quantum enhancement 
factor. The exemplary data correspond to the case of phase esti- 
mation using N photons in a Mach-Zehnder interferometer with 
5% losses in both arms. 

A or A' with probabilities p and 1 — p. The channels 
that cannot be decomposed into a convex combination 
of different channels (e.g. the unitary transformations) 
are called extremal. Note that while all interior points 
of the space of quantum channels are non-extremal, the 
boundary contains both the extremal as well as some 
non-extremal channels. 

We say that the family A v is classically simulated 38 if 
each channel is written as a classical mixture of the form 

K\P\ =Jdxp lp (x) A x [p], (4) 

where the unknown parameter enters only through the 
probability distribution p v of a random variable X that 
indicates which channel to pick from the set {A x } (see 
Fig. lb). If X\, . . . , Xn are independent hidden random 
variables used to generate the parallel channels we can 
rewrite the estimation problem as 

Since conditionally on the values of Xi the output state 
does not carry any information about <p, the estimation 
precision in the above scenario is at most equal to that of 
the classical problem of estimating <f> given N indepen- 
dent samples from the distribution p^ 

Hence, using the classical Cramer-Rao bound (CRB) 43 
we obtain a lower bound on the uncertainty of the original 
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Figure 3 | Local classical simulation. Schematic representation 
of a local classical simulation of a channel that lies inside the 
convex set of quantum channels. 



problem (1) 



Aip N > 



dx 



[d v P v {x)Y 



(5) 



where F c i is the classical Fisher information; see Methods 
for an alternative derivation. 

Provided that p v (x) satisfies regularity conditions re- 
quired in the derivation of CRB 43 , the SS of precision 
follows immediately. Moreover, any particular classical 
simulation yields an explicit SS bound on precision. If 
many CSs are possible, we obtain the tightest SS bound 
for A^jv by choosing the "worst" decomposition yielding 
minimal -F c i[7v], as shown below. On the other hand, HS 
is possible only when the above conditions are not sat- 
isfied. This happens in the decoherence-free case where 
A v are unitary channels which are extremal points of 
the SQC and the only admissible p v in (4) is the Dirac 
delta distribution being zero on all channels except A^. 
This yields F c \ =oo, makes the bound trivial and leaves 
the possibility of HS, which is consistent with results of 
decoherence-free metrology. 

The QFI at a given ipo depends only on the output 
state and its first derivative at <pq. This implies that any 
family of channels A v which "locally coincides" with the 
original one, i.e. 
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achieves the same maximum QFI in (3). It is therefore 
enough to consider the local classical simulations, i.e. any 
mixtures reproducing the channel and its first derivative 
at given tpo . As proven in Supplementary Methods, Sec- 
tion I, the CS with the smallest F c i can be constructed 
using two channels, {A + , A_}, that lie on the tangent line 
to the "channel trajectory" at the two outermost points 
situated on the boundary of the set of channels, sec Fig 3. 
The local CS around 4>o reads explicitly 



K[P\ =P% A +[p] +P V A_[p], 



(0) 



with 



A±[p)=A Vo [p] ±e±d v A v [p]\ pt 



e T ±{(p- tp ) 

£+ + £- 



Making use of Eq. (5) applied to the binary probability 
distribution , we obtain 



Aip N > 



N 



(7) 



To calculate the above bound it suffices to find the "dis- 
tances" e± of the channel from the boundary measured 
along the tangent line. For extremal channels e± = and 
the bound is not useful. For non-extremal channels, the 
above construction will yield a finite F c \ provided that 
e± > 0, i.e. A Va can be decomposed into a mixture of 
channels lying on the tangent. Channels which have this 
additional property will be called tp -non- extremal. They 
obey the standard precision scaling, and include all full 
rank channels (i.e. channels lying in the interior of the 
SQC) 38,39 . An explicit method for calculating e± and 
hence the bound for a general quantum channel is de- 
scribed in the Methods section. 



CHANNEL EXTENSION 

Even though the CS method is very general, there are 
interesting examples of (p-extremal decoherence models 
for which CS does not apply. In this case one can resort 
to the more powerful but less intuitive CE method. 

The action of the quantum channel A v can be de- 
scribed via its Kraus representation 44 : 

A lp [p]=Y<Ki{ip)pKt{ip), 

i 

with Kraus operators satisfying Yli-^-ii'P) ^iif) = 1- Al- 
though this representation is not unique, different sets of 
linearly independent Kraus operators are related by uni- 
tary transformations 



where u(<^>) is a unitary matrix depending on tp. 

An equivalent definition of the QFI has 
proposed 37 ' 39 



Fq[A, 



min F Q (\* V )), 



(8) 



been 



(9) 



where the minimization is performed over all ip differen- 
tiable purifications of the output state 

A>] = Tr s {|^X^|}- 

For pure input state, different purifications correspond to 
different equivalent Kraus representations of the channel, 
as in Eq. (8). For many quantum metrological models 37 
one can make an educated guess of a purification and de- 
rive excellent input independent analytical bounds pro- 
viding the correct asymptotic scaling of precision. Never- 
theless, the method may be cumbersome especially when 
the channel description involves many Kraus operators. 
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A simpler bound can be derived by exploiting the in- 
tuitive observation that allowing the channel to act in a 
trivial way on an extended space, can only improve the 
precision of estimation, i.e. 

maxFo [A v [p\] <m&xF Q [A v <g> I [p ext ]} ■ 

P Pcxt 

This leads to an upper bound on J-jy, which goes around 
the input state optimization, yielding 39 : 

T N < 4 min {N \\ ait \\ + TV (TV - 1) ||/J /? || 2 } (10) 

where ||-|| denotes the operator norm, the minimization 
is performed over all equivalent Kraus representations of 
A v , and 

«*=E^i(v)^i(v). Pk = '^Zfo>)Kii!P)> ( n ) 

i i 

where Ki{tp) = d v Ki(tp). For any given tp = tpo, equa- 
tion (10) involves only Ki and its first derivatives at tp . 
Moreover, the bound is insensitive to changing the Kraus 
representations with a tp independent u. Therefore, it is 
enough to parameterize equivalent Kraus representations 
in Eq. (8) with a hermitian matrix h which is generator of 
u(tp) = exp[— ih(tp — fa)]- This reduces the optimization 
problem (10) to a minimization over h. 

Since the SS of precision holds when Tn scales linearly 
with TV, the bound (10) implies that a sufficient condition 
for SS is to find h for which = 0, or equivalently 39 : 

i,j 1 

Here we go a step further and show that in this case we 
can obtain a quantitative SS bound 

Fn < 4TVmin||a!^|| , (13) 

where the minimisation runs over h satisfying (12), and 
can be formulated as a semi-definite program, as de- 
scribed in the Methods section. 

Moreover, it turns out that the bound resulting from 
finding the global minimum in Eq. (13) is at least as 
tight as the one derived using the CS method based on 
Eq. (5) (see Supplementary Methods, Section II) proving 
superiority of the CE over the CS method. 

EXAMPLES 

All examples of channels presented below are of the 
form A v [p] = A.[U v pU^\ , i.e. a concatenation of a uni- 
tary rotation encoding the estimate parameter tp and 
an 93-independent decoherence process. Consequently, 
Ki(ip) = KJJ V , where Ki(ip), Ki are the Kraus opera- 
tors of K v and A respectively. The most relevant models 
in quantum enhanced metrology belong to this class, but 



the methods presented may be applied to more general 
models as well. 

In what follows, we adopt the standard notation where 
1 is the 2x2 identity matrix and {<Tj} i _ 1 are the Pauli 
operators. We focus on two-level probe systems (qubits) 
sensing a phase shift modeled using a unitary U(tp) — 
cxp[(i<T3<^)/2] — rotation of the Bloch ball around the 
z axis. In the case of atomic clocks' frequency calibra- 
tion ip = 5u> ■ t with 6lu being the detuning between the 
frequency of the atomic transition and the frequency of 
driving field, while t the time of evolution. Even though 
in practice the parameter to be estimated is 8uj, we will 
consider that to be tp, in order to have a unified nota- 
tion for atomic and optical models. In the case of a two 
mode optical interferometer, U(ip) is the operator acting 
on a single photon state, accounting for the accumulated 
relative phase shift tp between the two arms of the inter- 
ferometer. 

We apply the methods to four decoherence processes 
encountered in quantum enhanced metrology: two-level 
atom dephasing, depolarisation, spontaneous emission 
and the photon loss inside the interferometer. These ex- 
amples will provide us with the full picture of the appli- 
cability of the methods discussed in the paper, as their 
cover all distinct cases from the point of view of the ge- 
ometry of SQC. 

Depolarization. Two level atom depolarization de- 
scribes an isotropic loss of coherence and may be visual- 
ized by a uniform Bloch ball shrinking, see Fig. 4a, where 
< t) < 1 is the final Bloch ball radius. Its description 
involves four Kraus operators 

*-=^*- {*-vtM (i4 » 

1 ) i=i. ..3 

which makes it an example of a channel lying in the 
interior of the SQC. Using CS method we infer the 
SS of precision, and calculate the "distances" from the 
boundary of the SQC e± = - r))(l + 3rj)/(2r]) (see 
Methods) resulting in the bound given in left column 
of Tab. Ia. Applying the CE method it is possible to 
further improve the bound and the result of the semi- 
definite optimization is shown in the right column of 
Tab. Ia. The only non zero elements of the correspond- 
ing optimal h are h 3 = h 30 = a/(1 — 77) (1 + 3r])/c and 
hi2 = -h 2 i = -i(l + 17) /c where c = 2(1 - 77) (1 + 2rj). 
To our best knowledge, the bound has not been derived 
before. 

Dephasing. Dephasing is a decoherence model of 
two-level atoms subject to fluctuating external mag- 
netic/laser fields. In graphical representation, it corre- 
sponds to shrinking of the Bloch ball in x, y directions 
with z direction intact, see Fig. 4b. The canonical Kraus 
operators read 44 

*° = ^¥«. *. = ^3, (15) 

where < r\ < 1 is the dephasing parameter. As it involves 
only two Kraus operators, it is not a full-rank channel and 




Figure 4 | Decoherence models. Graphical representation of decoherence models discussed in the paper. Two-level atom decoherence 
processes are illustrated with a corresponding shrinking of the Bloch ball, a, Depolarisation b, Dephasing c, Spontaneous emission 
d) Lossy interferometer with 77 being the power transmission coefficient. 



channel considered 



classical 
simulation (7) 



channel 
extension (13) 



depolarisation 



b dephasing 



spontaneous 



lossy 
interferometer 



(l-77)(l + 3»7) 



4^ 



l+»7-2i7^ 
2^P 



N/A 



N/A 



Table I | Precision bounds of the most relevant models in quan- 
tum enhanced metrology derived with two methods discussed 
in the paper. All the bounds are of the form A^at > , 
where constant factors are given in the table. Classical simula- 
tion method does not provide bounds for spontaneous emission 
and lossy interferometer, since these channels are ip-extremal. 
For the dephasing model it surprisingly yields an equally tight 
bound as the more powerful channel extension method. 



with < i] < 1. Interestingly, for all 77 this channel is 
extremal 45 , which means that the CS is not applicable. 
Nevertheless, the CE method of Eq. (13), can be easily 
employed. 

Substituting Ki(ip) into Eq.(12) we find that the gen- 
erator h is fixed to h = j^xz^Tj (°» — vtyi Consequently, 
there is no need for minimization over h and from 
Eq. (13) we automatically obtain an SS bound listed in 
Tab. Ic which to our knowledge has not been reported in 
the literature before. 

Lossy interferometer. In order to model loss in an 
optical interferometer, a third orthogonal state at the 
output, — vacuum — resulting from a loss of a photon, 
needs to be introduced. The decoherence channel on a 
single probe (single photon) is a map from a two to a 
three dimensional system: 




VT^T] 



K 2 = 



^ 




\/T=rj 



(17) 



lies on the boundary of the SQC. It is, however, a non- 
extremal and more importantly non-<^-extremal channel 
allowing for the CS construction with e± = \Jl — rf /r) 
(see Methods), yielding the bound given in Tab. Ib. 

Most importantly, the above bound, is exactly the 
same as the one derived by Escher et al. 37 or by the CE 
method, where the minimum in Eq. (13) corresponds to 
h=<7i/(2-\/l — T) 2 ). This proves that despite its simplic- 
ity the CS method may sometimes lead to bounds that 
are equally tight as the ones derived with much more ad- 
vanced methods even for channels lying on the boundary 
of the SQC. 

Spontaneous emission. The well known two-level 
atom spontaneous emission model is described by the 
Kraus operators 



K n = 



1 

y/fj 



K x = 



VT 








(16) 



where r\ is the power transmission coefficient for the light 
traveling through each of the two arms. Although the 
corresponding channel A v is non-extremal, it is unfor- 
tunately c/j-extremal and CS cannot be used. Still, CE 
method can be easily applied. The optimal bound (see 
Tab. Id) corresponds to h with non-zero elements h 0u = 
— hn — 2(1—^) ■ This bound is asymptotically equally 
tight to the best bounds known in the literature 35-37 
proving again that CE method despite its simplicity is 
able to provide powerful results in a straightforward man- 
ner. 

In optical interfcromctric applications it is common to 
use states of light with an unbounded number of photons 
such as coherent or squeezed states. At a first glance it 
is not obvious that the model considered in the paper 
covers these situations. Formally speaking the bound we 
have derived applies to input states with a total number 
of photons fixed to N. Notice, however, that in every op- 
tical experiment what is measured in the end are photon 
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numbers. If all phase reference beams are taken into ac- 
count we can regard quantum states of light as incoherent 
mixtures of states occupying different total photon num- 
ber sectors 46 ' 47 . From the point of view of metrology the 
QFI is then bounded from above by the weighted sum 
of QFIs for each of these sectors 30 . Within each sector 
we can easily apply our bound and since the bound is 
linear in N the effective bound will simply correspond to 
replacing N with N — mean number of photons in all 
relevant beams used in the experiment — and as such 
can be automatically applied to experiments involving 
coherent and squeezed light. 

In order to demonstrate the practical relevance of the 
bound it is instructive to compare its predictions with 
the actual quantum enhancement observed recently in 
the GEO600 gravitational wave detector 17 . A detailed 
theoretical analysis of the setup from the perspective of 
the derived bounds is underway. Still, by reducing the 
essential features of the setup to a simple Mach-Zchnder 
interferometer, we can give a preliminary estimate on 
how far is the actual experiment from the optimal per- 
formance. For the reported overall optical transmission 
77 = 0.62 the theoretically predicted maximal quantum 
enhancement amounts to a s/1 — rj = 0.62 factor reduc- 
tion in estimation uncertainty compared with the classi- 
cal l/y/r]N limit. The reported experimentally observed 
reduction was a factor of 0.67 which is an indication that 
the experiment operates close to the fundamental quan- 
tum limit and any significant improvement is possible 
only if the optical loss is further reduced. 



SUMMARY 

Assessing the impact of decohcrcnce on the maximum 
possible quantum enhancement is a crucial element in 
developing quantum techniques for metrological applica- 
tions. The tools developed in the paper allow for a direct 
calculation of bounds on the precision enhancement for 
arbitrary parameter estimation model where the deco- 
herence process acts independently on each of the probes 
and may be represented with a finite number of Kraus 
operators. While the CE method is more powerful and 
for the most relevant metrological models yields in the 
asymptotic limit of large number of probes the tightest 
bounds known in the literature, the CS method may fail 
to provide equally tight bounds but provides an intuitive 
geometric insight into the absence of asymptotic HS in 
the presence of decohcrcnce. From the derived bounds it 
is clear that if the HS was to be preserved for large num- 
ber of probes N the level of decoherence would have to de- 
crease with increasing N roughly as (1 — 77) rj 1/N . This 
gives an estimate on the regime in which the quantum en- 
hancement is quadratic as compared to the regime of con- 
stant factor improvement. This is clearly seen in Fig. 3 
where this transition appears around N ~ 1/(1 — rf) =20 
for r) = 0.95. Since in most metrology applications N is 
larger by several order of magnitude, we expect that the 



SS scaling provides a reliable bound for the optimal esti- 
mation precision. 

An important question that was not addressed in the 
paper is the saturability of the bounds. Clearly the CS 
method does not provide tight bounds in general as we 
have observed e.g. for depolarization channel. As to the 
CE method we are not able to prove that the bounds we 
have derived are saturable in the asymptotic limit. We 
should stress however, that these are lower bounds on 
estimation uncertainties and an upper bound can always 
be found by choosing a particular estimation method. 
Therefore if a particular strategy performs close to the 
derived bound, we can certify that it is near the fun- 
damental quantum limit, as illustrated by the results of 
GEO600 experiment . 

Note that even though the minimization over purifica- 
tion method 37 yields in principle a tight bound, in prac- 
tice there is no effective algorithm to find a global mini- 
mum and the only way to convince oneself that the one 
has hit the global minimum is again to show that the the- 
oretical lower bound on estimation uncertainty coincides 
with the performance of a particular estimation strategy. 

The methods discussed were focused on deriving use- 
ful bounds in the asymptotic regime of large number of 
probes. Still, the bounds are valid (though weaker) for 
any value N. We leave it for a future work to improve the 
bounds for finite N which seems to be possible by using 
the CE method and relaxing the /3j> = constraint. 

METHODS 

Proof of the classical simulation bound. In order to 
simplify the reasoning, let us focus on the classical simulations 
that are constructed using discrete sets of quantum channels, 
{A;}. Then, the considered channel's action of Eq. (4) can 
be rewritten as a (^-independent map acting on a larger input 
space 38 

A^[p]=^iW A 4p] = $b®0V]: (18) 

i 

where a v = Y^, p Vt i \et){ei\ represents a diagonal state in a 
basis, in which $ is defined via <&[g] = '^2, i A.i®Ei [g] with 
Ei[a] = (e«| <J |ej). In order to prove (5), we write the QFI for 
the N parallel use of the channel and bound it from above, 
i.e. 

F Q [Ar[p)]=F Q [$® N [ P ®*r]} 

<F Q [p9>a^] =F Q [a® N ] 

= NF Q {<T v \=NF cl ]p v \, (19) 

exploiting the monotonicity of the QFI under any parameter 
independent quantum map 48 , here $® JV . 

Calculation of the classical simulation bound. The 

geometry of the space of channels and more specifically the ip- 
extremality are best viewed by using the Choi-Jamiolkowski 
isomorphism 49 ' 50 . Given a quantum channel A : £.(Hin) i-> 
C(H \it) acting from the space of density matrices on "Hi n to 
density matrices on H ou t, one defines Pa = A <g> 
where |$) =X^ii Win i s a maximally entangled state 
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in Hin ® Hin- A is a physical channel (i.e. trace pre- 
serving, completely positive map) iff Pa is a positive semi- 
definite operator, satisfying Tr-H out {-PA} = 1. If {Ki}i are 
the Kraus operators of the A channel, we can write explicitly 
Pa = Ei \Ki)(Ki\, where \K t ) =K t ® 

We can now say that, the channel A v is not (^-extremal, if it 
is possible to find a non-zero e, for which P\ v ± eO^Pa^ > 0. 
See Supplementary Methods, Section III for an alternative 
formulation of the <£>-extremality condition and its relation 
to the well known non-extremality condition due to Choi 50 . 
In practice, if we want to make most out of the bound in 
Eq. (7), we need to find the maximum values of e±, for which 
Pa v ±£± d v PA v are still positive semi-definite operators. This 
is a simple eigenvalue problem and therefore the bound can 
be obtained immediately. 

Taking the dephasing model as an example, the Choi- 
Jamiolkowski isomorphism of the corresponding A v channel 
Pa v =Ei=i \KiU v }(KiU v \ has a simple form: 



Pa, = 



1 77 \ 





1 



(20) 



It is easy to check that Pa v + e d v P\ v > provided |e| < 
y/l — V 2 /v> hence using Eq. (7) we arrive at the bound given 
in Table I. 

Channel extension method as a semi-definite pro- 
gram. Here we show that the minimization problem in 
Eq. (13) can be formulated as a simple semi-definite program. 
Let the channel K v be a map from a di to a d-z dimensional 
Hilbert spaces with Kraus representation involving k linear 
independent Kraus operators (cfe X di matrices). Consider 
the following block matrix: 



A = 



( V~tt dl 
ki 



V~tt d2 
o 



V Kk-i 





Vtt d2 
o 



(21) 



Vitd 2 J 



where 1^ is a d x d identity matrix. Positive semi-definiteness 
of the matrix A is equivalent to the condition 



*j i = Y,kli( i <ti dl . 



Minimizing the operator norm \\a^\\ is thus equivalent to min- 
imizing t subject to A > 0. Taking into account Eq. (12) the 
problem takes the form: 

min t, subject to: A > 0, ^ hijKjKj = i ^ k\ K q . (22) 



Since Ki = Ki — ihij Kj the matrix A is linear in h and 
the problem is thus a semi-definite program with the resulting 
minimal t being the minimal operator norm Uo^H- For the 
purpose of this paper we have implemented the program using 
the CVX package for Matlab 51 . 
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Appendix A: Optimal tangent simulation 

We prove that the optimality of CS is reached by choos- 
ing two channels, which are as far as possible from A Vo on 
the tangent, but still lie in the set of physical maps. Con- 
sider a CS being a mixture composed of channels that lie 
on the tangent line to the channel trajectory. Pick a chan- 
nel Ao situated a position xo, Ao = A Vo + xq d v K v \ , 
which contributes to the mixture with p° . Let us now 
remove this channel from the mixture and replace it by 
two channels A± at positions x±, X- < xq < x+, with 
mixing probabilities . For the total probability to sum 
up to one and the new mixture to be locally equivalent 
to the previous one, the following quantities as well as 
their first derivatives over tp should equalize at ifo: 



P^+P V =P° V , P~?pX+ +P V x- =f>° ab- 
solving the above set of equations we get 



(Al) 



± _ o x ± 

P<Po ~ Pvo _ , 



X 



dipP v | 



dyp'lWo 



x± - x 



x± 



(A2) 

We now want to show that the described operation can 
only decrease the Fisher information F c \. For a moment 
consider a situation in which channels A± were not used 
in the original mixture and only appeared at the expense 
of the removed Aq channel. The change of F c \ is given by 



5F cA = 



Pv 



Pv 



(A3) 



ip=ip 



and by Eq. (A2) it must be zero. 

If, on the other hand, channels A± were already in 
use in the original mixture, we can artificially separate 
their initial part from the one originated at the expense 
of taking away probability from the point xq- This for- 
mally corresponds to introducing additional distinguish- 
ing symbol that would split the instances of the classical 
random variable x± into two groups, depending on their 
origin. By the same reasoning as before, the F c \ must 
remain constant during probabilities redistribution. Fi- 
nally, removing the distinguishing symbol we can only re- 
duce the F c i, what proves that indeed replacing A with 
A± in a local classical simulation can only decrease the 
F c i. 

This reasoning may now be applied to an arbitrary 
probability distribution of channels lying on the tangent 
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line. We can only decrease the F c i by deleting one of the 
channels and redistributing its probability to the ones 
lying at the line ends. Repeating this procedure for all 
channels on the tangent line we conclude that the optimal 
(yielding minimum F c {) simulation is the one composed 
of two channels situated at the intersection of the tangent 
line with the boundary of the set of all quantum channels. 

The missing point in the above reasoning is to prove 
that it is not possible to decrease F c \ further by in- 
cluding channels that lie outside the tangent line. We 
have proved this fact for the class of channels A<p[p] = 
A(U v pU^) — concatenation of a unitary parameter en- 
coding and a decoherence map. The sketch of the proof 
is the following. Notice that the values of parameters Xj 
appearing in Eq. (A2) will now be obtained by perform- 
ing a projection of the channel onto the tangent line. The 
only way, in which the use of this larger class of channels 
could lead to lower F c i, would be the situation when after 
projecting a channel onto the tangent line, we find our- 
selves outside the physical region of all quantum maps. 
The effective xi lies then further apart than when consid- 
ering only channels on the tangent line. This seems like 
a possible way to increase the product e + e_. However, 
taking one channel, e.g. A+, to the "top" of the tangent 
line requires taking the complementary channel, i.e. A_, 
to the "bottom". Since the decoherence structure is tp- 
invariant and the set of channels is convex, the gain in 
the £ + will be compensated by at least an equal loss in 
the distance e_, which will not make the product e+e— 
larger, and therefore will not decrease the F c \. 



Moreover, notice that 



K I 



(B3) 




m? , [d P -] 2 



(B4) 



Substituting the above inequality to Eq. (13) we see that 
Fn < NF C \, where F c \ is just the Fisher information for 
the local CS. This shows that the bound derived from CS 
will never be tighter than the one derived from the CE 
method optimized over all Kraus decompositions. 



Appendix C: tp non-extremality 

We prove below that the condition on <p non- 
cxtrcmality, which requires the existence of a non-zero 
epsilon such that 



P hv ±ed^P Av >0, 



(CI) 



is equivalent to the statement that there exist a non zero 
Hermitian matrix p.ij such that 



d v PA v =E^« \ K *)( K o 



(C2) 



Appendix B: Channel extension vs. classical 
simulation method 



If a channel can be locally simulated at tpo, than up to 
the first order in tp — tpo we have 

A v [p] = p+ A+[p] + P - A_[p] 



E^9^?' 



where 



{*,} 




structed Kraus operators satisfy Eq.(12), since 



The con- 

incc 

% = iE k l R i = 1 ( E K K PKf+ (Bl) 

q ^ i 

+ E*v*7 + *r) = 1 +Pv) 1 = °< ( B2 ) 

i ' 
where we have used V y : p+ + p~ = 1 . 



Assume that (CI) holds and recall that = 
\Ki)(Ki\. Since the operator on the l.h.s. is positive 



E 



Ki){Ki\ ± e J2 \Ki){Ki\ + \KitKi 



> 



(C3) 

for any vector \tp). In order to prove (C2), it is enough to 
show that \Kj) can be written as linear combinations of 
K t . If this was not the case for one of vectors, e.g. \Kj), 
we would additionally need a vector \Lj) that is orthog- 
onal to the the space spanned by then taking 
IV') = \/p\Kj) + exp(i£)\/l — p\Lj) positivity condition 
(C3) leads to 

pimK^ 2 ± y/p{\-p)e (e^miKj) + cc) > (C4) 

For any nonzero s and nonzero (L^\K*) we can always 
find some ^ and p small enough, so that the l.h.s. is 
negative. This leads to a contradiction, hence (C2) must 
hold. 

For the opposite direction, assume (C2) holds and sub- 
stitute the formula for dP^ v into the l.h.s. of (C2). We 
need to show that there exist e such that 



A± = E \ K *x K i i ± e E \ k m k 3 i ^ °- 



(C5) 
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Defining matrices 



1 ± e/Lt, note that, for e 



small enough, they are positive semi-definite. Hence, 
we can take their square root and construct \K^) = 



\Kj). Then, A ± = EJ^X^I, which is 



clearly positive semi-definite. ■ 

As a last remark, observe that tracing out (C2) over 
Hout space yields 



which is the Choi condition for non-extremality 45 . This 
reflects the trivial fact that if a channel is tp non-extremal, 
then it must not be extremal. 



(C6) 
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